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ABSTRACT 

With the aim of studying magnetic effects in time- distance helioseismology, 
we use the first-order Born approximation to compute the scattering of acoustic 
plane waves by a magnetic cylinder embedded in a uniform medium. We show, 
by comparison with the exact solution, that the travel-time shifts computed in 
the Born approximation are everywhere valid to first order in the ratio of the 
magnetic to the gas pressures. We also show that, for arbitrary magnetic field 
strength, the Born approximation is not valid in the limit where the radius of 
the magnetic cylinder tends to zero. 

Subject headings: MHD, Magnetic Fields, Scattering, Waves, Sun: Helioseismol- 
ogy, Sun: Magnetic Fields 

1. Introduction 



Time-distance helioseismology (jDuvall et al.lll9931 ) has been used to measure wave travel 



times in and around magnetic active regions and sunspots to estimate sub surface flows and 



wave-speed perturbations (e.g. iDuvall et al.l Il996l : iKosovichev et al.ll2000l ). A challenging 



problem is to estimate the subsurface magnetic field from travel times. In order to do so, 
one must understand the dependence of the travel times on the magnetic field. 



As discussed by e.g. ICallyl ( 120051 ) the interaction of acoustic waves with sunspot magnetic 
fields is strong in the near surface layers. As a result, the effect of the magnetic field on the 
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travel times is not expected to be small near the surface. Deeper inside the Sun, however, 
the ratio of the magnetic pressure to the gas pressure becomes small, and it is tempting 
to treat the effects of the magnetic field on the waves using perturbation theory. The 
hope is to eventually develop a linear inversion to estimate the subsurface magnetic field 
from travel times measured between surface locations that are free of magnetic field. Of 
particular interest is the search for a magnetic field at the bottom of the c onvec tion zone. 
Such a linear inversion scheme has been proposed by iKosovichev fc Duvalll (119971 ) for time- 
distance helioseismology using the ray approximation, but it needs to be extended to finite 
wavelengths. 

As a first step, in this paper, we consider the scattering of small amplitude acoustic 
plane waves by a magnetic cylinder embedded in a uniform medi um. This simp le problem 
has a known exact solution for arbitrary magnetic field strengths (1 Wilson! 1 1 98 Ol ) . The first- 
order Born and Rytov approximations have proved useful in the context of time-distance 
helio seismoloKy to model the effects of small local perturbations in sound speed and flows 
(e.g. iBirch et al.ll200ll ; I Jensen fc Pijpera 120031 : iBirch fc Felderl I2004J ) . Here we use the Born 
approximation to compute the scattering of a wave by a weak magnetic field. The validity 
of the Born approximation is not a priori obvious in this case, since the magnetic field allows 
additional wave modes. Because we have an exact solution, however, we can study the 
validity of the linearization of travel times on the square of the magnetic field. We note that 
the problem of the scattering of waves by a non-magnetic cylinder with a sound speed that 
differs from the surrounding medium was investigated by iFan et al.l (119951 ) . 



The outline of the paper is as follows. In Section [2] we specify the problem and write the 
equations of motion for small amplitude waves. In Section [3] we review the exact solution 
to the scattering problem. In Section H] we apply the first Born approximation to obtain 
the complex scattering amplitudes. In Section we show that the Born approximation is 
an asymptote of the exact solution in the limit of infinitesimal magnetic field strength. In 
Section [6] we compare travel times computed exactly, in the Born approximation, and in the 
ray approximation. In Section [7] we provide a brief summary of our results and also discuss 
the limit when the magnetic tube radius tends to zero. 
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2. The Problem 

2.1. Governing equations 

We start with the ideal equations of magnetohydrodynamics. The equations of conti- 
nuity, momentum, magnetic induction, and Gauss' law for the magnetic field are: 

0, (1) 

0, (2) 
0, (3) 
0, (4) 

where D t = d t + v ■ V is the material derivative, p the density, v the velocity, p the pressure, 
and B the magnetic field. For the sake of simplicity, we use the simple energy equation 

pC v D t T + pV -v = 0, (5) 

where T is the temperature, C v = R/(j — l) the uniform specific heat at constant volume, R 
the gas constant, and 7 the ratio of specific heats. This equation neglects all forms of heat 
losses. In addition we use the ideal gas equation of state 

p = pRT. (6) 



Dtp + pV-v = 

pD t v + Vp- ^-(VxB)xB = 

d t B -Vx(vxB) = 

V B = 



2.2. Steady background state 

We consider a magnetic cylinder with radius R and uniform magnetic field strength B t 
embedded in an infinite, otherwise uniform, gravity free medium with constant density po, 
gas pressure po, and temperature T . We use a cylindrical coordinate system (r, 6, z) where 
r is the radial coordinate, 9 is the azimuthal angle, and z is the vertical coordinate in the 
direction of the cylinder axis. We denote the corresponding unit vectors by f, 0, and z. All 
steady physical quantities are denoted with an overbar. In particular, we have 

B = B t Q(R - r)z, (7) 
p = p t Q(R-r)+p e(r-R), (8) 
p = p t Q(R-r)+p e(r-R), (9) 

where the Heaviside step function is defined by 6(r) = if r < and 6(r) = 1 if r > 0. 
The density and pressure inside the tube are p t and p t respectively. We assume that there 
is no mean flow in this problem, i.e. v = 0. 
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We choose to study the case where the background temperature is the same inside and 
outside the magnetized region, i.e. T = T t = T . As a result, the sound speed, c = (7-RT) 1 / 2 , 
is constant everywhere. This choice is motivated by our desire to restrict ourselves, as much 
as possible, to the study of the effect of the Lorentz force on waves, rather than the effect of 
a sound speed variation. Pressure balance across the magnetic tube boundary implies 

pt + B?/8n = po, (10) 

where p t is the background gas pressure inside the tube. The density inside the tube is given 
by p t = PoPt/Po, as the temperature is the same inside and outside the tube. 



2.3. Linear waves 

We want to study the propagation of linear waves on the steady background state defined 
above. Toward this end, we expand each physical quantity that appears in equations ([I])-© 
into a time- varying component, denoted by a prime, and the steady component, denoted 
with an overbar. For example, we write p = p + p'. After subtraction of the steady state, we 
obtain: 

d tP ' = -V •(/*/), (11) 
pdtv' + Vp' = -J- \(VxB')xB + (VxB)xB'] , (12) 

47T 

d t B' = Vx (v'xB) , (13) 
V B' = 0. (14) 

The linearized energy equation, in combination with equation (TTTT) and the linearized equa- 
tion of state, may be simplified to 

d t p' - c 2 <V = ^— *W • Vp, (15) 

7 

which describes adiabatic wave motion. 

As we study linear waves on a steady background, we can consider one temporal Fourier 
mode at a time. The magnetic field B and all other background quantities do not depend on 
z. Thus, a wave with a z dependence of the form e lkzZ will have the same z dependence after 
interacting with the magnetic cylinder. As a result, we study solutions where the pressure 
fluctuations are of the form 



p'(r, z, t) = p{r) exp(ik z z — iut), 



(16) 
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where r = (r, 9) is a position vector perpendicular to the tube axis. All the other wave 
variables, p', v', and B' are written in the same form as equation f|T6|) . Quantities with a 
tilde only depend on r. 



Exact Solution 



For the sake of completeness, we briefly review an exact solution obtained by IWilson 



(Il980l ) to equations ffTTj) - ffT5j) . We consider a plane wave incident on the magnetic tube, with 
pressure fluctuations of the form 



PincO) = Pexp(ik ■ r) 



(17) 



where P is an amplitude and k is the component of the wave vector perpendicular to the 
tube axis. In order for the incident wave to be a solution to the non-magnetic problem, the 
horizontal wavenumber k = \\k\\ must satisfy 



k{u) = v^Vc 2 - k 2 . 



18) 



In the rest of this paper, unless otherwise stated, we will use k to denote k{uj). In cylindrical 
coordinates, this pla ne wave can b e expanded as a sum over azimuthal components (index 
m) according to (e.g. lBogdanlll989l ): 



Ph 



P J2 i m J m (kr)e m 



(19) 



m=— oo 



where J m denotes the Bessel function of order m and is the angle between k and r. 



The total wave pressure, hydrodynamic plus magnetic, and the radial v elocity must be 
continuous across the tube boundary. Applying these boundary conditions, IWilsonl (119801 ) 
showed that the total pressure wave field is 



p(r) 



PEm imR mJm{ktr)e im ^ r<R 
fc + ^Em'^W^ r>R 



(20) 



where H m = Hm is the Hankel function of the first kind of order m. The quantity k t is the 
horizontal wavenumber inside the tube, given by 



2„2\ 



kia 



(1 + a 2 1 'c 2 )(lu 2 - s 2 k 2 z ) 



1/2 



(21) 
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1/2 

where a = B t / y/Airp t is the Alfven wave speed and s = ac (a 2 + c 2 ) is the tube velocity. 
The coefficients A m and B m are given in Appendix |A] This exact solution is valid for arbi- 
trarily large values of B t and R. We note that the gas pressure fluctuations are discontinuous 
at the tube boundary 

Not all azimuthal components m contribute at the same level. Figure [T] shows the 
amplitudes of the A m as a function of m for two different values of the tube radius, R. In 
this particular example, po — 5 x 10~ 7 cgs, c = 11 km/s, and B = 1 kG are fixed. For 
R = 0.5 Mm, which is less than the wavelength (kR = 0.86), only the m = 0, ±1 azimuthal 
components contribute. For tubes with larger radii the higher order components have larger 
amplitudes, as can be seen in the case R = 2 Mm (kR = 3.43). 



4. First Born approximation 

Magnetic effects cause perturbations to both the steady background state and the wave- 
field. In this section we use the Born approximation to derive an approximate solution 
to equations (fTTD - (fT5|) based on the assumption that magnetic effects are small (see e.g. 
Rosenthal 19951 ) . The Lorentz force is quadratic in the magnetic field. As a result we intro- 



duce a small parameter that is second order in the magnetic field. We choose to expand all 
quantities in powers of the small dimensionless parameter 

p2 

(22) 



In this expansion framework the magnetic field appears at order e 1//2 . In particular, we write 
the steady background magnetic field as 

B = e 1/2 B 1 . (23) 

The magnetic field causes a shift, epx, in the steady component of the density inside the 
tube relative to the steady component of the density outside the tube, p : 

p = p + e Pl . (24) 

Likewise, we write the steady component of pressure as 

p = p + ep 1 . (25) 

These changes are related to the magnetic field through equation (fTDl and the equation of 
state: 

Pi = -^0(R-r), (26) 
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T 




Azimuthal order, m 

Fig. 1. — Coefficients |v4 m | as a function of m for R = 2 Mm and = 0.5 Mm (thick lines). 
Only the m > values are shown since |A_ m | = \A m \. The magnetic field is B = 1 kG, 
a;/27r = 3 mHz, and k z = 0. The coefficients A m are negligible for m > 5 in the case 
R = 2 Mm and for m > 2 in the case i? = 0.5 Mm. The open circles are the Born 
approximation to these coefficients (|A^ orn |, see section 0J. 
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Pi = -^e(fl-r). (27) 

We expand each of the wave-field variables into an incident component (subscript "inc" ) 
and a scattered component (subscript "sc"): 

V = pL + wL (28) 

P' = Anc + ePso (29) 

v' = < c + e< c , (30) 

B' = e^B' sc . (31) 



By inserting the above expansions into equations fllljl - fTToTl and retaining the terms of order 
e we obtain 

- iwfb + Po V • v' sc = -V-(pi< c ), (32) 
-iup v' 8c + Vp' 8c = iup l v' hlc + —(VxB 1 )xB BC 

+i(VxBjxB 1) (33) 

-iu (p' sc - c 2 p' sc ) = v[ nc - Vpi, (34) 

-iuB' sc = Vx« c x J B 1 ). (35) 

The terms on the right-hand side of the above equations act as sources for the scattered waves: 
this is the Born approximation. Writing all wave variables in the form of equation (fl6|) and 
using the fact that the magnetic field is solenoidal, the above equations reduce to a forced 
Helmholtz equation for the (k z ,uj) Fourier component of the scattered pressure field, p sc : 

(A r + k 2 )p sc (r) = S(r), (36) 

where Ar is the two-dimensional Laplacian with respect to r and the source function S(r) 
is given by 

S(r) = 1 ^p L S(r - R)d r p inc (r) 

( ~ k2 -{A r -kl)[Q{r-R)p mc {r)} 



LU 2 
C 2 



2u< 



(A r - 3k 2 ) [8{r - R)d r p- mc (r)] . (37) 



The first term in 5* is due to the density jump at the tube boundary and the other two terms 
are due to the direct effect of the Lorentz force on the wave. For the incoming wave, pi nc , 
we take the same plane wave as in the exact solution (Eq. |17j). 
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Fig. 2. — a) Real and b) imaginary parts of the scattered pressure field in the Born ap- 
proximation (dashed line) and the exact solution (solid line) . In this case the incoming wave 
is of the form p inc = Jo(kr) (m = 0) and we used B = 1 kG, R = 2 Mm, k z = 0, and 
u)/2n = 3 mHz. 
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The solution to the inhomogeneous Helmholtz equation fl36l) is 

P«(r) = J J G(r\r')S(r')dr' (38) 

where G(r\r') is the Green's function defined by 

(A r + k 2 ) G(r\r') = 6(r - r') (39) 
and explicitly given by (e.g. iMorse fc Ingardlll986l ) 

oo 

G ( r IO = ~\ E ^ m (A:r > )J m (fcr < )e^- e '), (40) 

m=— oo 

where r = (r, 9), r' = (r', 0'), and r> = max(r, r') and r< = min(r, r'). 

We insert this expression for G into equation (138]) and use integration by parts as 
appropriate. The solution can be written in terms of integrals over bilinear combinations of 
Bessel and Hankel functions. These integrals can be evaluated using equations fIBlj) and flB2j) . 
Upon simplification, we find that the pressure of the scattered wave is 



ep sc (r) = P ^2 



f C m J m {kr) - t^J'Jkr) r<R 

\ Af r rH m (kr) r>R, 1 1 

where the coefficients A^ rn and C m are given in Appendix O 

Figure [2] shows the real and imaginary parts of the Born approximation for the scattered 
pressure field, ep sc , for the case pi nc = Jo(kr), i.e. when m = and P = 1. In this example, 
e = 0.13. On the same figure, we also show the exact calculation of the scattered pressure 
field obtained by subtracting p inc from the complete exact pressure field, p, computed as in 
Section O The amplitude of the scattered wave is of order e of the amplitude of the incoming 
wave, as expected. The Born approximation is everywhere accurate to about 10%. 

For a direct comparison with the exact solution (Eq. [20]). we remind the reader that 
the total pressure wave field in the Born approximation is given by p = p mc + ep sc , according 
to equation 



5. Born tends to the exact solution as e — > 



In this section we show that to first order in e, the exact solution (Sec. [3]) and the Born 
solution (Sec. H]) are identical. We expand the exact solution (Eqs.[20], |Alj . and |A2j ) in a 
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Taylor series up to first order in e. To do this, we use 



a 2 



c 2 ~ 1 - 7 e/2 (42) 

which, together with equation (l2Tj) . gives the first-order perturbation to the wavenumber 
inside the tube, 



A; 4 (e) = A;(l-e/2) + 0(e 2 ). (43) 

Let us first consider the exact scattering coefficient A m outside the tube (r > R) given by 
equation flAlj) . Denoting the numerator and denominator of A m by iV and D respectively 
and performing a Taylor expansion, we obtain 



N = i e 



2 1,2 

(7 + 2^)J' m {kR)J m (kR) + kRJ 2 m {kR) 
us 1 



2 

-kRJ m ^(kR)J m+1 (kR)\ + 0(e 2 ) (44) 



and 



£> = J m (kR)H' m (kR) - J' m {kR)H m {kR) + 0(e) 
7r/cit 

Hence, the exact and Born coefficients outside the tube match to first order in e: 

A m = Al°™ + 0(e 2 ). (46) 

Similarly, it can be demonstrated that the coefficient B m that gives the exact total wave 
field inside the tube (r < R) is 

B m -l = C m + 0(e 2 ). (47) 

The minus one on the left-hand side comes from the fact that the B m coefficient relates to 
the full wavefield, whereas C m is for the scattered wavefield only. Together, equations (|4"6"1) 
and (j4"7j) imply that the Born approximation is identical to the exact solution outside and 
inside the magnetic tube, to first order in e. Figure [3] shows the fractional error r] m = 
l^m° m — An|/|An| as a function of e for m = 0, m = 1, and m = 2. We see that the 
fractional error of the Born approximation tends to zero as e tends to zero. 



6. Travel Times 

In this section, we study the interaction of solar-like wave packets with the magnetic 
cylinder. The aim is to compare seismic travel-time shifts computed in the Born approxima- 
tion and from the exact solution. For the sake of simplicity, we fix k z = 0. Using Cartesian 
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0.0 0.1 0.2 

Perturbation strength, e 

Fig. 3. — The fractional error r\ m = \A®° rn — A m \/\A m \ as a function of e for m = (solid 
line), m = 1 (dotted line), and m = 2 (dashed line) for the case R = 2 Mm, o;/27r = 3 mHz, 
and /c z = 0. 
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coordinates r = (x, y) for the horizontal plane, we choose an incoming Gaussian wavepacket 
propagating in the +x horizontal direction: 

POO 

[ nc (r, t) = / e -("-^) 2 /(^ 2 ) CO s[k(u)x - tut] du. (48) 



Pi: 

where u^/2n = 3 mHz is the dominant frequency of solar oscillations and <r/27r = 1 mHz is 
the dispersion. Since we chose k z = 0, we have k{uS) — lu/c. The wavepacket is centered on 
the magnetic tube at time t = 0. The scattered wave packet can be calculated, exactly or in 
the Born approximation, from the previous sections. 

The three panels in Figure H] show snapshots of the incoming and scattered pressure fields 
at time t = 8.9 min. The parameters of the steady background at infinity are po — 5x 10 -7 cgs 
and c = 11 km/s, which are roughly the conditions at a depth of 250 km below the solar 
photosphere. The incident wavepacket is shown in Figure Hk- Figure Hb shows the scattered 
wave which results from a 1-kG magnetic flux tube with radius R = 0.2 Mm. The smaller 
tube produces relatively more back scattering than the large tube in comparison with the 
forward scattering. The amplitude of the scattered wave is roughly three orders of magnitude 
smaller than the incoming wave. Figure Hb shows the scattered wave for a larger tube radius 
of 2 Mm; the scattering is dominantly in the forward direction and has an amplitude only 
an order of magnitude smaller than the incident wave. 

We now define the travel-time shifts that are caused by the magnetic cylinder. By 
definition, the travel-time shift at location r is the time 5t(r) which minimizes the function 

X(t) = J dt' \p'(r, t') - p[ nc (r, H - t)f , (49) 

where p' is the full wavefield that includes both the incident wavepacket and the scattered 
wave packet caused by the magnetic field. The travel-time shifts can be computed in this 
way for either the exact solution or the Born-approximation. 

In addition, it is also interesting to c ompa re with the ray approximation as given by 



equation (14) from iKosovichev fc Duvalll (119971 ). In our case, where k ■ B = and the 



magnetic field strength is constant inside the tube, the ray approximation becomes 5t(r) = 
—L(r)a 2 /(2c 3 ) where L(r) is the path length through the tube along the ray which goes 
from coordinates (— oo,y) to r = (x,y). 

Figure [5] shows travel-time shifts resulting from a flux tube of radius 2 Mm and field 
strength of 1 kG (e = 0.13). Figure shows the exact, Born-, and ray- approximation 
travel times as a function of x at fixed y. Inside the flux tube, both the Born- and ray- 
approximation travel times reproduce the exact travel times at a good level of accuracy. As 
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x increases to the right of the tube, wavefront healing (e.g. iNolet fc Dahlenl 120001 ) is seen 
in the exact and Born approximation travel times. Wavefront healing, however, is not seen 
in the ray approximation travel times. Figure \5jp shows the travel times as a function of y 
at fixed x = 10 Mm. The Born approximation reproduces the exact travel times to within 
20%. The ray approximation does not capture finite wavelength effects and does not capture 
the basic behavior of the travel times; it can be inaccurate by many orders of magnitude 
for kR <C 1. We note that, in Figure [5j the contribution of the density jump (first term in 
Eq. [27]) to the travel-time shifts is negligible compared to the contribution from the Lorentz 
force. 



Discussion 



We have computed, in the first Born approximation, the scattering of acoustic waves 
from a magnetic cylinder embedded in a homogeneous background medium. We showed 
that in the limit of weak magnetic field, the Born approximation to the scattered wavefield 
is correct to first order in the parameter e = B 2 /4npc 2 . For typical values of the solar 
magnetic flux, the Born approximation should be good at depths larger than a few hundred 
km below the photosphere. The condition e << 1 is satisfied for a 1-kG magnetic fibril at a 
depth of 250 km (e ~ 0.1) and for a 10 5 G magnetic flux tube at the base of the convection 
zone (e ~ 10~ 7 ). Since the err ors introduced by the Rytov and Born approximations are very 



similar (e.g. IWoodwardl Il989l ) . we suspect that a travel-time shift computed in the Rytov 



approximation would also tend to the exact solution as e tends to zero. 



Near the photosp 



rere, e is not small. It has been suggested by many authors (e.g. 



Lindsey fc Braunl l2004f) th a t in t his case the Born approximation will fail. An exception 
is the claim by iRosenthal Jl995h that the Born approximation will remain valid for kG 
magnetic fibrils in the limit where the radius of the magnetic element is much smaller than 
the wavelength. We wish to test this last statement in our simple problem. 



Assuming k z 
for all e we have 



for the sake of simplicity and taking the limit kR 
A* orn { l + fl-^e ifm = 0, 



lim 

kR^O A r 



l + (l- 7 /2)e 
1 - 7e/4 



otherwise. 



0, we find that 



(50) 



This shows that the Born approximation is not valid in the limit of small tube radius. Figure[6] 
shows the ratio A^ Tn / A m , for < m < 5, as a function of R when e = 1 and k = 3.7 Mm -1 . 
We see that, in the limit of small kR, the fractional error in the Born approximation is of 
order e (the absolute error is of order e 2 ). The Born approximation applied to completely 
evacuated solar magnetic fibrils in the photosphere is likely to be invalid by roughly a factor 
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of two. Note that the sign of the relative error in A^ Tn is different for m = and m > 0. 
The sensitivity of travel times to local perturbations in internal solar pr operties can be 



descri bed through linear sensitivity functions, also called travel-time kernels. iGizon fc Birch 



(120021 ) gave a general recipe for computing such travel-time kernels using the B orn approx 



i matio n, which has been applied to the case of sound-speed perturbations by iBirch et al. 



(120041 ) . The present work suggests that travel-time kernels for the subsurface magnetic field 
will be useful for probing depths greater than a few hundred km beneath the photosphere, 
at least in the case when the travel times are measured between surface points that are not 
in magnetic regions. One should be careful, however, not to draw definitive conclusions from 
the simple model we have studied, given the complexity of the real solar problem. 

A.C.B. was supported by NASA contract NNH04CC05C and thanks the Max-Planck- 
Institut fur Sonnensystemforschung for its hospitality. 
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A. Exact solution coefficients 



The coefficients A m and B m are 

- ( l + iS) k iJ' m (hR)J m {kR) + (l - Jm{k t R)J' m (kR) 



and 
B,, 



2i 



(! + If) ^J' m {hR)H m {kR) 



nkR 



k 
kt 



1 + 



7 a 
2^ 



1 - 



A- 2 



a 2 fc? 



J m {k t R)H' m {kR) 



(Al) 



J' m (k t R)H m (kR) - ., 



a 2 k 2 



J m (k t R)H' m (kR) 



(A2) 

In equations (lAlj) and (1A2|) . the functions J' m and H' m denote the first derivative of J m and 
H m = Hm respectively. 



B. Useful Integrals 



In order to compute scattering amplitudes in the Born approximation, we used (IWatson 



1944 chap. 5) 



(Bl) 



This preprint was prepared with the AAS IATgX macros v5.2. 
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and 

„2 



j x'H m (kx')J m (kx') dx' = — [2J m (kx)H m (kx) - J m ^ 1 (kx)H m+1 (kx) - J m+ i(kx)H m _ l (kx)} 

(B2) 



C. Born approximation coefncients 

The coefficients A® orn and C rn for the Born solution are 



a Born _ inkR 

■ft-™ — c 



m 



2 t,2 

7 + ■>'— ) J' m (kR)J m (kR) + kRJ 2 m {kR) - />•/?■/,„- . ( A //)./„,. ,j i /, /?) 



(CI) 



and 



^ m = - C ^--e^ 7 + 2^) J'JkR)H m (kR) 



UJ 2 4 \' LO 2 

\2 



-e \ ' [2J m (kR)H m (kR) - J m ^(kR)H m+1 (kR) - J m+1 (A; J R)iJ m _ 1 (A; J R)j;C2) 
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Fig. 4. — Plots of the pressure field of the unperturbed (a) and scattered (b - c) wavepackets 
at time t = 8.9 min after the unperturbed wavepacket has crossed the center of the magnetic 
cylinder. The wavepacket parameters are described in the text. The wavevector is normal 
to the axis of the magnetic cylinder and in the +x direction. The strength of the magnetic 
perturbation is e = 0.13. In panel (b) the tube radius is R = 0.2 Mm and in panel (c) is the 
tube radius is R = 2 Mm. The circles outline the cross-section of the tube. Notice that the 
gray scales are different in each panel. The backscattered wave is more prominent for the 
tube with the smaller radius (panel b). 
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Fig. 5. — Local travel-time shifts St(r) caused by the magnetic cylinder (e = 0.13). The 
travel times are measured at positions r in a plane perpendicular to both the cylinder axis. 
The incoming wavepacket, which moves in the +x direction, is the same as in figure Hh- The 
radius of the tube is R = 2 Mm and the tube axis is (x, y) = (0, 0), as shown in figure Hfc. 
In both panels the heavy solid line is the exact travel-time shifts, the circles are the Born 
travel-time shifts, and the light line gives the ray approximation. The left panel shows the 
travel-time shifts as a function of x at fixed y = 0. The right panel shows the travel-time 
shifts as a function of y at fixed x — 10 Mm. The Born approximation is reasonable for 
this value of e. The ray-appr oximation does not cap ture finite-wavelength effects and fails 
to describe wavefront healing (INolet & Dahlenll200 



r 
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Fig. 6. — Ratio A^ rn /A m in the complex plane at fixed e = 1 and k z = 0. The ratio is plotted 
for varying values of the tube radius in the cases m = (thick line), m = 1 (thin line), and 
2 < m < 5 (dashed lines). The big circles show the limit kR given by equation ( l50l) . If 
the Born approximation were correct for small tube radii, the big circles would coincide with 
the cross. The small and medium-size circles are for kR = 1 and kR = 1/2 respectively. 



